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The new three year WMAP data seem to confirm the presence of non-standard iarge scale features 
in the cosmic microwave anisotropy power spectrum. While these features may hint at uncorrected 
experimental systematics, it is also possible to generate, in a cosmological way, oscillations on large 
angular scales by introducing a sharp step in the inflaton potential. Using current cosmological data, 
we derive constraints on the position, magnitude and gradient of a possible step. We show that a 
step in the inflaton potential, while strongly constrained by current data, is still allowed and may 
provide an interesting explanation to the currently measured deviations from the standard feature- 
less spectrum. Moreover, we show that inflationary oscillations in the primordial power spectrum 
can significantly bias parameter estimates from standard ruler methods involving measurements of 
baryon oscillations. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The recent three year results from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite 0, H H H 
have further confirmed with an extraordinary precision 
the inflationary paradigm of structure formation in which 
primordial fluctuations are created from quantum fluctu- 
ations during an early period of superluminal expansion 
of the universe |5j, |y, [7| • 

Indeed, soon after the WMAP data release, a number 
of authors investigated the possibility to discriminate be- 
tween several single-field inflationary models u sing this 
new, high quality, dataset @, 0, US, EH E3, [13, E^l 
One of the main conclusions of these works is that some 
inflationary models, such as quartic chaotic models of the 
form V((f>) ~ \<fi 4 , may be considered ruled out by the 
current data, while others, such as chaotic inflation with 
a quadratic potential V{<j>) ~ m 2 <j) 2 , are consistent with 
all data sets. 

One important assumption in these analyses (apart 
for [l^) is that the inflaton's potential is featureless, 
i.e., there is no preferred scale during inflation and the 
primordial power spectrum of density perturbations in 
Fourier fc-space can be well approximated by a power 
law k n , where the spectral index n is almost scale inde- 
pendent. The main prediction of these models is that the 
anisotropy angular power spectrum should be "smooth" 
and not show features in addition to those provided 
by the baryon-photon plasma oscillations at decoupling 
within the framework of the standard ACDM model of 
structure formation. 



The current WMAP data is in very good agreement 
with this hypothesis: several non-standard features in 
the anisotropy angular power spectrum detected in the 
first year data have now disappeared thanks to the longer 
integration time of the observations and better control of 
systematics (see [l|). 

However, features in the large scale anisotropy spec- 
trum are still present in the new release. Moreover, 
some of the cosmological parameters derived from the 
new WMAP data, like, for instance, the low value of 
the variance of fluctuations erg, appear in tension with 
those derived by complementary data sets. It is therefore 
timely to investigate a larger set of inflationary models 
and to consider a cosmological origin of these unexpected 
features. 

A departure from power law behavior of the primordial 
power spectrum could be caused by a chan ge o f the initial 
conditions, due to trans-planckian physics [16|,[l7| or un ~ 
usual initial field dynamics [l8L Il9| or by some brief vio- 
lation of the slow roll conditions during inflation [23, ■ 
We will investigate a model of the second type where 
features in the temperature and density power spectra 
arise due to a step-like change in the potential parame- 
ters, as proposed by J^]. A sharp step in the inflaton 
mass, caused, e.g., by a symmetry breaking phase tran- 
sition, generates indeed fc-dependent oscillations in the 
spectrum of primordial density perturbations. 

The goal of our paper is to make use of the recent three 
year WMAP data (WMAP3) and other datasets to con- 
strain the possibility of a step feature in the inflaton po- 
tential. For this purpose we adopt the phenomenological 
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model proposed by Adams et al. [22J, where a step fea- 
ture is added to the chaotic inflationary potential in the 
following way: 

V((f>) = im 2 2 ^l + ctanhT^p^ , (1) 

where m is an overall normalization factor, c determines 
the height of the step, d its gradient and b is the field 
value on which the step is centered. Previous phe- 
nomenological studies of the same j23| or other oscilla- 
tory features |2J, |25| have been limited to the first year 
WMAP data, and in general a full analysis varying also 
all cosmological parameters is still missing and is the ma- 
jor result of this work. 

The paper is organized as follows: In Sec. II we briefly 
review step-inflation models. In Sec. Ill we describe our 
analysis method. In Sec. IV we present our results and, 
finally, in Sec. V we derive our conclusions. 

II. INFLATION MODELS WITH A STEP IN 
THE POTENTIAL. 

Inflationary models with a step can naturally arise in 
theories with many interacting scalar fields, e.g., in super- 
gravity models. In general, these models contain several 
flat directions in field space and thus offer the possibility 
to have multiple inflationary phases separated by phase 
transitions |26 |. or even inflation with a curved trajectory 
in field space |27l I2H , |2^ | . In the last case the presence 
of additional active scalar degrees of freedom generates 
not only the adiabatic mode of curvature perturbations, 
but also the isocurvature one. Since the data do not 
seem to require an isocurvature component [3(j, we will 
restrict ourselves to the case where the phase transition 
does not appreciably change the rolling direction corre- 
sponding to the inflaton field and the energy density is 
always dominated by a single field. Also, we will investi- 
gate the simplest scenario and assume that the sole effect 
of the phase transition is to change the parameters of the 
Lagrangian for the inflaton field, in particular its mass. 

Consider a hybrid inflationary potential of the type 

V = V + \rnl4? + ^ - A/ 2 ) 2 + A 2 0V , (2) 

where <fi is the inflaton field, while ip is a hybrid field 
that takes a vacuum expectation value during inflation 
when the inflaton reaches the critical value 4>1 = I* 1 
the usual case of hybrid inflation, this transition is so 
strong that it stops the inflationary phase. But, if the 
coupling A is sufficiently small, the back-reaction of ip is 
too weak and inflation continues. On the other hand, the 
parameters of the Lagrangian change their value and for 
example the inflaton effective mass becomes 

m 2 ff (0)=m 2 + A 2 (V, 2 )(0). (3) 

In this scenario, even if the classical inflaton is nearly 
unperturbed, the inflaton perturbations are affected and 



the primordial power spectrum is modified. 

It is well known that a step in the mass generates os- 
cillations in the primordial spectrum and this can be 
described analytically in the WKB approximation [3l| . 
The behavior of the inflaton mass is determined by the 
dynamics of the phase transition and the growth of the 
hybrid field fluctuations, which become tachyonic after 
the critical point. This growth depends on the classi- 
cal inflaton field motion, but is in general so fast that 
ip reaches the minimum in a very small number of e- 
foldings [22, IH, |3 Ull ■ The inflaton mass is therefore 
reasonably well approximated by a hyperbolic tangent, 
and so we take 

m 2 ff (» ~ to 2 ^1 + ctanh (J^f^J ) • ( 4 ) 

Here, we see that the parameter m is an average infla- 
ton mass and b is of the order of the critical value <p c . 
The other two parameters determine the duration of the 
transition and the strength of the effect on the inflaton's 
mass. Note that we work in Planck units, so all dimen- 
sional quantities like b and d should be multiplied by Mp 
in order to obtain their value in physical units. In this 
article we restrict ourselves to the case of chaotic infla- 
tion, where the mass term determines both the classical 
dynamics of the inflaton and the behavior of the per- 
turbations. A more general discussion will be left for 
a longer publication j3g. We therefore assume that Vq 
and also its change due to the phase transition are com- 
pletely negligible. We also discuss here only positive c 
values; negative c is also allowed, but is restricted to be 
very small to avoid the presence of another minimum in 
the potential away from = 0. 

Since we cannot rely on the slow roll approximation 
for a generic choice of parameters, we integrate the equa- 
tions for the background and for the modes numerically 
as discussed in detail in p2| . The equations for the infla- 
ton field and the Hubble parameter in Planck units are 
simply 

4> + 3H4>+V'(<t>) = (5) 

3ff 2 = y + n0)- (6) 

We assume slow roll as the initial condition for </> 3> <fi c 
and solve the evolution numerically until the end of in- 
flation in order to determine the number of e-foldings 
between b and the end of inflation. 

The equation for the Fourier components of u — —zlZ, 
the curvature perturbation [37j . takes the usual form 

< + (k 2 - ^) «* = (7) 

where z — a<j>/H is given by the background dynamics 
and the primes and dots denote derivatives with respect 
to conformal time and physical time, respectively. Using 
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FIG. 1: Effects of a step in the potential on the power spec- 
trum of curvature perturbations. Here we show the primor- 
dial spectrum for the two best fit points, corresponding to 
b = 14.81, c = 0.0018, d = 0.022 (dashed line) and b = 14.34, 
c = 0.00039, d = 0.006 (solid line). 



the equations of motion for the classical field, we have in 
general 



— = 2a z H z 
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where V' and V" are the derivatives of the potential with 
respect to the inflaton field. In the step model, this quan- 
tity can deviate substantially from the slow-roll expecta- 
tion z" I z ~ 2a 2 H 2 at the time of the transition. Wc 
solve equations © to numerically using a Bulirsch- 
Stoer algorithm for free field initial conditions at an ini- 
tial time when k 2 3> z"/z is satisfied. 

Once we know the solution for the mode k, we can 
determine the primordial power spectrum 
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evaluated when the mode is well outside the horizon. Our 
results are stable with respect to changes in the exact 
time when we set the initial conditions and when we com- 
pute the spectrum. 

The resulting spectra as a function of k for different 
parameters are shown in Figure ^ Essentially, the spec- 
trum shows a power-law behavior with a superimposed 
oscillation. 

How will the four parameters of our model affect the 
shape of the spectrum? The overall normalization of Viz. 
is proportional to m 2 , b determines the wavelength at 
which the feature appears and the maximum amplitude 
of the oscillations is roughly proportional to c. Gener- 
ally, the dominant contribution to z" / z comes from the 
V" term and is proportional to c/d 2 , so the range of k af- 
fected by the feature depends on the square root of c/d 2 . 
Note that away from the step, the slow roll conditions are 
satisfied and the spectrum recovers the usual power law 
form with spectral index given by n s = 1 — 2/N ~ 0.96 
for a number of e-foldings N equal to 50, which is what 



we assume for our analysis. The value of the spectral 
index is the same before and after the step since it does 
not depend on m 2 in m 2 </> 2 models. Also, for values of 
the parameters where the slow-roll conditions are always 
satisfied (i.e., small values of c/d 2 ), the spectrum does 
not show a full oscillation, but a dip at the scales cor- 
responding to the transition. So even in this case it is 
not so well approximated by the usual power-law with a 
constant spectral index. 



III. CMB ANALYSIS 

We compare the theoretical model described in the pre- 
vious sections with a set of current cosmological data 
by making use of a modified version of the publicly 
available Markov Chain Monte Carlo (MCMC) package 
cosmomc |3Sj . 

We sample an eight-dimensional set of parameters. 
Four of them determine the primordial power spectrum, 
namely the b, c and d parameters of the step-inflation 
model as described in the previous section and the over- 
all normalization of the primordial power spectrum As 
(equivalent to m 2 as discussed earlier). The remaining 
four cosmological parameters are the physical baryon and 
CDM densities, io h = Q, b h 2 and oj c = fl c h 2 , the ratio of 
the sound horizon to the angular diameter distance at 
decoupling, S and finally, the optical depth to reioniza- 
tion, r. Furthermore, we consider purely adiabatic initial 
conditions, impose flatness and neglect neutrino masses. 

We include the three year data yj (temperature and 
polarization) using the likelihood routine for supplied 
by the WMAP team and available at the LAMBDA web 
site. 1 We marginalize over the amplitude of the Sunyaev- 
Zel'dovich signal. The MCMC convergence diagnostics 
are done on four chains using the Gelman and Rubin 
"variance of chain means" / "mean of chain variances" R 
statistics for each parameter, demanding that R—l < 0.1. 
Our 2D constraints are obtained after marginalization 
over the remaining "nuisance" parameters, again using 
the programs included in the cosmomc package. 

In addition to CMB data, we also consider the con- 
straints on the real-space power spectrum of galaxies 
from the Sloan Digital Sky Survey (SDSS) 

We restrict the analysis to a range of scales over which 
the fluctuations are assumed to be in the linear regime 
(k < 0.2 /i/Mpc). When combining the matter power 
spectrum with CMB data, we marginalize over an ad- 
ditional nuisance parameter the dark versus lumi- 
nous matter bias. Furthermore, we make use of the 
HST key project measurement of the Hubble parameter 
H = 100ft. km s^Mpc" 1 H2 b y multiplying the like- 
lihood by a Gaussian centered around h = 0.72 with a 
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FIG. 2: Marginalised likelihood (solid line) and projection of 
the likelihood distribution (dotted line) for the b parameter 
in the case of WMAP only. Two peaks for b at b = 14.3 and 
b — 14.8 are clearly visible. The one at b = 14.8 provides a 
good fit to the low £ WMAP glitches. It is evident that the 
likelihood function is far from gaussian in this direction. The 
difference between the two curves is caused by a volume effect 
when integrating over the other parameter directions. 



standard deviation a = 0.08. Finally, we impose a top- 
hat prior on the age of the universe, 10 < to < 20 Gyrs, 
and a Gaussian prior on £lf,/i 2 centered around 0.022 
with a standard deviation of 0.002 from BBN constraints, 

cf. Ref. jai. 

We demand that the feature appear at a wavelength 
to which our data is sensitive, so our analysis will 
be limited to the interval 13.5 < b < 16. Apart from 
that we also impose logarithmic priors on the other 
step parameters: log c £ [— 6, — 1], logd e [— 2.5, — 0.5] 
and logc/d 2 G [—5,3] and flat priors on the cosmologi- 
cal parameters. 

As it turns out, the likelihood distribution C has a 
rather odd shape and some of the interesting features 
are at low likelihoods. In order to improve mixing and 
get a better coverage of the low likelihood regions, we 
sample C 1 ^ 3 instead of C (i.e., we use "heated" chains at 
T = 3). 

As a measure of the performance of the step model, we 
compare its best fit x 2 with the best fit x 2 of a reference 
model, which we take to be the "vanilla" 6 parameter 
(ilbh 2 , Q c h 2 , 6 S , r, As and n s ) power law ACDM model. 



IV. RESULTS 

Let us first consider the WMAP dataset alone. In Fig- 
ure we plot the mean likelihood for the b parameter 
which determines the position (scale) of the step in the 
potential. If some value of b is preferred by the data, 
then it would hint at the presence of a feature. As we 
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FIG. 3: This plot shows the temperature anisotropy angular 
power spectrum of the best fit step model (WMAP only, solid 
line) and, for reference, the best fit 6 parameter power law 
ACDM model (dashed line) . The dotted line shows the effect 
of a feature near b — 14.8 for WMAP data only, i.e., the 
"local" best fit at the lower peak in Figure [5] 



can see, the mean likelihood distribution clearly indicates 
two maxima for the b parameter at b = 14.3 and b — 14.8, 
respectively. The feature at b = 14.8 is able to produce 
a good fit to the WMAP low-£ glitches (yielding an im- 
provement of Ax 2 — 5 over the vanilla model) and agrees 
with the results of Ref. [23 for fixed cosmological pa- 
rameters and WMAP first year data. Interestingly, the 
authors of J^j find an oscillating feature at roughly the 
same scale by reconstructing the primordial power spec- 
trum from the first year data. The minimum at b = 14.3 
provides oscillations on smaller scales (A^ 2 — 7), beyond 
the second peak in the anisotropy power spectrum, see 
Figure 

It is interesting to project the likelihood function onto 
the (6, logc) plane (Figure 0] upper panel). First of all, 
we see that for a range of values (14.3 < b < 15.5) a 
region of step models with logc < —3 is ruled out at 
99% confidence level. This disfavored region corresponds 
to the region in k space that is better sampled by the 
WMAP data, and where, therefore, the data provide 
the strongest constraints. Secondly, the two aforemen- 
tioned maxima in 6-space can again be seen in the 2D 
projection. These maxima are close to the two bound- 
aries (large and small scales) of the region sampled by 
WMAP and centered around amplitude logc ~ —3. We 
find that the WMAP polarization and cross temperature- 
polarization data are rather insensitive to the presence of 
features. 

The d parameter is not well constrained by the data 
due to a degeneracy with c. Instead, we consider the 
constraints in the (b, logc/d 2 ) parameter space (Figure^] 
bottom panel). This parameter is as well constrained as 
c and, again, the presence of two maxima for b is evident. 
Also, the maxima are at values of c/d 2 of order 1, where 
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FIG. 4: Mean likelihood and marginalized likelihood con- 
tours in the (6, c) and (6, c/d 2 ) planes at 8% and 99% c.l. for 
WMAP data only. The peaks comprise less than ten per cent 
of the total volume of the likelihood function. 



the slow-roll conditions are strongly violated; values of 
c/d 2 <C 0.1 correspond, on the other hand, to the usual 
slow-roll m 2 (j) 2 inflation and cannot be excluded by the 
data. 

Note that in the (6, c, <i)-subspace of parameter space 
the likelihood function is very oddly shaped and can by 
no means be approximated by a multivariate Gaussian. 
As a consequence, the likelihood at the boundaries of this 
subspace is generally not negligible (see, e.g., Figure Ell- 
in fact, the likelihood function will have a large plateau 
of constant likelihood in regions where the step model 
cannot be distinguished from the m 2 cf) 2 chaotic inflation 
model, either because the step is too small or too smooth 
(small c, large d) or because the feature appears at wave- 
lengths the data is not sensitive to (b too small or too 
large). Apart from the two peaks, we also find a valley 
that can be excluded at a high confidence level. How- 
ever, since the plateau may contain a significant fraction 
of the total volume of the likelihood distribution, con- 
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FIG. 5: Likelihood contours in the (b, c) and (b, c/d 2 ) planes 
at 68% and 99% confidence level for CMB+SDSS. 

straints derived from marginalization will be dependent 
on the priors on 6, c and d. 

For our choice of priors and WMAP data only, we find 
that the peak regions contribute about 8% of the total 
volume, so, from a Bayesian standpoint, the WMAP data 
alone do not require the presence of a feature. 

FigureElindicates that the best fit model has a feature 
at a range of wavelengths where the WMAP data are 
limited by large systematic errors. It is therefore inter- 
esting to enquire whether the inclusion of other data sets 
which are more sensitive at small scales will corroborate 
this result. To this end we add small scale CMB data 
from the ACBAR, BOOMERANG, CBI, MAXIMA and 
VSA experiments [IS H USEE! and the SDSS large 
scale structure data. 

These data sets probe mainly smaller scales and are 
not sensitive to the large scale feature at b = 14.8. 

We find that including these small scale data improves 
the constraints on the oscillations quite significantly: 
A\ 2 — 15, to which the SDSS data alone contribute 
about 6. As we can see from Figure (bottom panel), 
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FIG. 6: Top: temperature anisotropy angular power spec- 
trum with small scale CMB data. 
Bottom: galaxy power spectrum and the SDSS data. 
The solid lines depict the best fit step model, the dashed line 
the best fit reference model. 



step-induced oscillations in the matter power spectrum 
seem to provide a much better fit to the SDSS data. 

In Figure El we plot the likelihood contours in the 
(6, logc) plane at 68% and 99% c.l. for the CMB+SDSS 
case. Adding the SDSS data increases the statistical sig- 
nificance of the maximum, the peak near (b ~ 14.3 , 
logc ~ —3.5, logc/d 2 ~ 0) now contains about 70% of 
the total volume of the likelihood function, the likelihood 
at the boundaries of parameter space is suppressed and 
hence, the results are much less prior dependent. Also, 
we can rule out a much larger chunk of parameter space 
for b < 14.5 due to the increased sensitivity of the data 
on the corresponding scales. 

The variance test convergence stats using last half 
chains are 0.0007 for b, 0.0038 for logc and 0.0036 for 
logc/d 2 , showing a robust convergence of the chains. 

We also considered data from the 2dF Galaxy Redshift 
Survey [U and we found that the CMB+SDSS results 
are stable under the inclusion of the 2dF dataset. 

Of course, oscillatory behavior in the observed data 



could well have its origin in uncorrected or unidentified 
systematic effects such as a scale dependent bias. How- 
ever, as indicated by our results, the presence of multi- 
ple step-like features in the inflaton potential (expected 
from, e.g., supergravity or M-theory models |4!j) is also 
a viable solution. 



It is important to check if the step in the potential sug- 
gested by the CMB+SDSS analysis has some impact on 
the estimation of the remaining cosmological parameters. 
We find no correlations between b, c, d and the cosmolog- 
ical parameters As,0 s and r. This is clearly due to the 
fixed spectral index the model has away from the feature. 

On the other hand some correlation is present between 
c and the baryon and cold dark matter energy densi- 
ties r2fc/i 2 and £l c h 2 . Figure shows the confidence level 
likelihood contours in the (fib/i 2 , c) and (il c h 2 , c) param- 
eter space. A deeper step in the potential (larger c) has 
the effect of making the data more compatible with a 
lower baryon density and a higher cold dark matter den- 
sity. While the effect is small, it is interesting to note 
that the baryon density derived from the WMAP data in 
the framework of the standard model is generally larger 
than that predicted by standard big bang nucleosynthe- 
sis and measurements of the primordial deuterium abun- 
dance (see, e.g., 
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FIG. 7: Correlations between c and the cosmological param- 
eters Q b h 2 and Q, c h 2 for CMB+SDSS. Plotted are the 50%, 
68% and 99% confidence levels. The best fit step model peak 
prefers a slighly lower value of Q, b h 2 than the plateau region 
where the feature is insignificant. One can also see that near 
the best fit region, the larger c the smaller the corresponding 
value of Qhh 2 . 



MIMICKING BARYONIC OSCILLATIONS VI. CONCLUSIONS 



Recently, a detection of oscillations in the correlation 
function of the Luminous Red Galaxies sample of the 
SDSS has been reported in If those oscillations 

are the imprint of primordial acoustic oscillations in the 
primeval baryon+photon plasma, then they may pro- 
vide a standard ruler at the redshift of the survey and 
a very powerful tool for testing the late time evolution 
of the universe and, ultimately, the appearance of dark 
energy. Clearly, if oscillations in the primordial spectrum 
are present they may mimic baryonic oscillations |5lj at 
different scales and drastically change the estimation of 
cosmological parameters. 

As a qualitative example, we plot in Figure |S| the cor- 
relation functions for a model with no baryons, the stan- 
dard ACDM model (Q& = 0.05) and a model with no 
baryons but with a step in the inflationary potential. As 
we can see, a model with no baryons but with oscilla- 
tions in the primordial spectra can reproduce the ob- 
served data very well. A more detailed and quantitative 
analysis will be presented in a future paper |36| . 
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FIG. 8: Correlation functions for a model with no baryons, 
the standard ACDM model (fit = 0.05) and a model with 
no baryons but with a step in the inflationary potential. The 
data points are from the LRGS sample of [ 



The new three year WMAP data seem to confirm the 
presence of non-standard large scale features on the cos- 
mic microwave anisotropy power spectrum. While these 
features may hint at uncorrected experimental systemat- 
ics, a possible cosmological way to generate large angular 
scale oscillations is to introduce a sharp step in the in- 
flaton potential. By making use of current cosmological 
data we derive constraints on the position, magnitude 
and gradient of a possible step in the inflaton potential. 
Our conclusion is that such a step, while strongly con- 
strained by current data, is still allowed and may provide 
an interesting explanation to the current measured devi- 
ations from the standard featureless spectrum at low I. 

Surprisingly though, the combination of all CMB data 
sets with the SDSS data seems to prefer a feature at 
small scales, which could mimic the effect of baryonic 
oscillations and reduces the best fit value of f2&. Note that 
for this it is sufficient to have a minute change, of order 
0.1%, in the inflaton mass parameter, but a relatively 
fast one. 

It is an open question if such a sharp step can be re- 
alized in realistic inflationary models and if the effect is 
physical. In general we can exclude the presence of strong 
features with c > 0.003 in the observable range. 

Future experiments like PLANCK will provide better 
measurements of the polarization and cross temperature- 
polarization spectra, providing an important check for 
possible non-standard features. 



Acknowledgments 

We would like to thank S. Sarkar for useful discussions. 
JH thanks Y. Wong for discussions and comments on 
the manuscript. LC and JH acknowledge the support 
of the "Impuls- und Vernetzungsfonds" of the Helmholtz 
Association, contract number VH-NG-006. 



[1] D.N. Spergel et ffl/.,|arXiv:astro-ph/0603449| 

[2] L. Page et al, |arXiv:astro-ph/0603450| 

[3] G. Hinshaw et al, arXiv:astro-ph/0603451 

[4] N. Jarosik et al, arXiv:astro-ph/0603452 

[5] AA. Starobinsky, JETP Lett. 30 (1979) 682 [Pisma Zh. 

Eksp. Teor. Fiz. 30 (1979) 719]. 
[6] VF. Mukhanov and G. V. Chibisov, JETP Lett. 33 

(1981) 532. 

[7] J.M. Bardeen, P. J. Steinhardt, and M. S. Turner, Phys. 
Rev. D 28 (1983) 679. 

[8] L. Alabidi and D.H. Lyth, |arXiv:ast ro-ph/0603539"1 

[9] H. Peiris and R. Easther, |arXiv:ast ro-ph/0603587 



R. Easth er and H. Peiris, |arXiv:astr o-ph/060 4'214] 
[10] A. Lewis, arXiv:astro-ph/0603753 

[11] U. Seljak, A. ST osar and P. McDonald, 

arXiv:astro-ph/0604335 
[12] J. Magueijo and R.D. Sorkin, arXiv:astro-ph/0604410 
[13] D. Parkinson, P. Mukherjee and A.R. Liddlc, 

|arXiv:astn>p h/0605003 C. Pahud, A.R. L iddle, 

P. Mukherjee and D. Parkinson, arXiv:astro-ph/0605004 
[14] W.H. Kinney, E.W. Kol b, A. Melchiorri and A. Riotto, 

arXiv:astro-ph/0605338 PRD in press. 
[15] J. Martin and C. Ringeval, arXiv:astro-ph/0605367 
[16] R.H. Brandenberger and J. Martin, Mod. Phys. Lett. 



8 



A 16 (2001) 999 [arXiv:astr o-ph/0005432 ; J. Martin 
and R. Brandenberger, Phys. Rev. D 68 (2003) 063513 
arXiv:hep-th/0305161 . 

R. Easther, B.R. Greene, W.H. Kinney and G. Shiu, 

Phys. Rev. D 66 (2002) 023518 |arXiv:he p-th/0204129 . 

CP. Burgess, J.M. Cline, F. Lemieux and R. Holman, 

JHEP 0302 (2003) 048 arXiv:hep-th/0210233 . 

C.R. Contaldi, M. Peloso, L. Kofman and A. Linde, 

JCAP 0307 (2003) 002 |arXiv:astro- ph/0303636 . 

L. Kofman, G.R. Blumenthal, H. Hodges and J.R. Pri- 

mack, ASP Conf. Ser. 15 (1991) 339. 

A. A. Starobinsky, JETP Lett. 55 (1992) 489 [Pisma Zh. 

Eksp. Teor. Fiz. 55 (1992) 477]. 

J. A. Adams, B. Cresswell and R. Easther, Phys. Rev. D 
64 (2001) 123514 arXiv:astro-ph/0102236 . 
H.V. Peiris et al, Astrophys. J. Suppl. 148 (2003) 213 
arXiv:astro-ph/0302225 . 

J. Martin and C. R ingeval, JCAP 0501 (2005) 007 
arXiv:hep-ph/0405249 . 

M. Kawasaki, F. Takahashi and T. Takahashi, Phys. Lett. 
B 605 (2005) 223 arXiv:astro-ph/0407631 . 
J. A. Adams, G.G. Ross and S. Sarkar, 
Phys. Lett. B391 (1997) 271 |arXiv:hep-p h/9608336 and 
Nucl. Phys. B503 (1997) 405 |arXiv:hep -ph/9704286]. 
V.F. Mukhanov and P. J. Steinhardt Phys. Lett. B422 
(1998) 52-60 arXiv:astro-ph/9710038 . 
A. A. Starobinsky, S. Tsujikawa and 
J. Yokoyama Nucl. P hys. B610 (2001) 383-410 
arXiv:astro-ph/0107555 . 

and B.J.W. van Tent 
19 (2002) 613-640 



5! Groot Nibbelink 
Class. Quan t. Gr ay. 

arXiv:hep-ph/0107272 . 
M. Beltra n et al, Phy s. 

arXiv:astro-ph/0501477 . 
P. Hunt and S. Sarkar, Phys. 

arXiv:astro-ph/0408138 . 
G.N. Felder, J. Garcia-Bellido, P.B. Greene, L. Kofman, 
A.D. Linde, I. Tkachev Phys. Rev. Lett. 87 (2001) 011601 



Rev. D 71 (2005) 063532 
Rev. D 70 (2004) 103518 



arXiv:hep-ph/0012142 . 
[33] T. Asaka, W. B uchmiiller and L. Co vi Phys. Lett. B510 

(2001) 271-276 |arXiv:hep-p h/0104037l. 
[34] E.J. Copeland, ST Pascoli and A. Rajantie 

Phys. Rev. D65 (2002) 103517 arXiv:hep-ph/0202031 . 
[35] J. Garcia-Bellido, M. Garcia Perez and 

A. Gonzalez-Arroyo P hys. Rev. D67 (2003) 103501 

arXiv:hep-ph/0208228 . 
[36] L. Covi, J. Hamann, A. Melchiorri, A. Slosar and I. Sor- 

bera, work in progress. 
[37] E.D Stewart and D. H. L yth, Phys. Lett. B 302 (1993) 

171 arXiv:gr-qc/9302019 . 
[38] A. Lewis and S. Bridle, Phys. Rev. D 66 (2002) 103511 

(Available from http://cosmologist.info). 
[39] M. Tegmark, A.J.S. Hamilton and Y.Z.S. Xu, 

Mon. Not. Roy. Ast ron. Soc. 335 (2002) 887 

|arXiv:astro-ph/01 1 1575| . 
[40] 3! Cole et al [The 2dFGRS Collaboration], 

Mon. Not. Roy. Ast ron. Soc. 362 (2005) 505 

[arXiv:astro-ph /0501 174| . 
[41] W.L. Freedman et al, Astrophys. J. 553 (2001) 47. 
[42] A. Shaf ieloo and T. S ouradeep, Phys. Rev. D 70 (2004) 

043523 [arXiv:astro-ph/0312174| . 
[43] J. H. G oldstein et of., Astrophys. J. 599 (2003) 773 

[arXiv:astro-ph/0212517]. 

[44] C.J. MacTavish et a<., |arXiv:a stro-ph/0507503] 

[45] A.C.S. Readhead et al, Astrophys. J. 609 (2004) 498 

[461 R- Stompor et al. [MAXIMA Collaboration], 

arXiv:astro-ph/0309409 
[47] C. Dickinson et al, Mon. Not. Roy. Astron. Soc. 353 

(2004) 732 [arXiv :astro-ph/0402498^ 

[48] B. Fields arid S. Sarkar, |arXiv:astro-ph/060i5l"4l 
[49] A. Ashoorioon and A. Krause, |arXiv:he p-th/0607001 
[50] D.J. Eisenstein et al, A strophys. J. 633 (2005) 560 

[arXiv:astro-ph/0501 171 1 . 
[51] X. Wang, B. Peng, M. Li, X.L. Chen and X. Zhang, Int. J . 

Mod. Phys. D 14 (2005) 1347 arXiv:astro-ph/0209242 . 



